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ABSTRACT: The gravitational Dirichlet problem - in which the induced metric is fixed 
on boundaries at finite distance from the bulk - is related to simple notions of UV cut¬ 
offs in gauge/gravity duality and appears in discussions relating the low-energy behavior 
of gravity to fluid dynamics. We study the Einstein-Maxwell version of this problem, in 
which the induced Maxwell potential on the wall is also fixed. For flat walls in otherwise- 
asymptotically-flat spacetimes, we identify a moduli space of Majumdar-Papapetrou-like 
static solutions parametrized by the location of an extreme black hole relative to the wall. 
Such solutions may be described as balancing gravitational repulsion from a negative-mass 
image-source against electrostatic attraction to an oppositely-signed image charge. Stan¬ 
dard techniques for handling divergences yield a moduli space metric with an eigenvalue 
that becomes negative near the wall, indicating a region of negative kinetic energy and 
suggesting that the Hamiltonian may be unbounded below. One may also surround the 
black hole with an additional (roughly spherical) Dirichlet wall to impose a regulator whose 
physics is more clear. Negative kinetic energies remain, though new terms do appear in the 
moduli space metric. The regulator-dependence indicates that the adiabatic approximation 
may be ill-defined for classical extreme black holes with Dirichlet walls. 
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1 Introduction 

Boundary conditions are an important part of any physical problem. For gravitational 
problems they are often imposed by restricting the falloff of the metric at infinity. Such 
asymptotic boundary conditions are appropriate for isolated physical systems. But there 
are also applications in which it is interesting to study the dynamics of a finite region of 
spacetime. 

A well known example is the study of thermodynamics for “gravity in a box” (see 
e.g. [1, 2]). In this regard it is interesting to note that simple black-hole-in-box calculations 
based on stationary black holes indicate that such regulated systems have bounded entropy 
and energy, as one would expect for held theories in finite volume with a UV cut-off, 
although the temperature diverges as the black hole horizon approaches the wall of the box. 
The apparent bound on area follows from the requirement that the black hole fit inside a 
box of fixed size, while the apparent bound on energy is the statement that Schwarzschild 
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Figure 1 . At the linearized level the Dirichlet wall will repel uncharged black holes, as can be seen 
using the method of images. However, if these black holes are extremal and oppositely charged, the 
electrostatic attraction precisely cancels the gravitational repulsion. 


black holes have finite energy in the limit where their size approaches that of the box [3]; 
see however section 5 of [4] for further comments. 

Taking the limit where the box walls approach the black hole horizon has been argued 
to lead to a controlled version of the black hole membrane paradigm [5, 6]. In particular, 
Bredberg et al. [7, 8] found in this limit that gravitational dynamics reduced to that of a 
fluid living on the cutoff surface (see also [9, 10]). 

In addition, the walls of our box can also be used to model the conhniug potential of 
AdS space. Indeed, in the AdS/CFT context it is common practice [11-16] to regulate AdS 
calculations by restricting to a finite volume V and fixing the metric on dV, thus imposing 
what from the CFT perspective is a surprising Lorentz-invariant UV regulator. See also 
[7, 17, 18] for related discussions of Wilsonian renormalization. 

The above discussions highlight the need for a more complete understanding of grav¬ 
itational dynamics in the presence of such boundaries. We take preliminary steps in this 
direction below by investigating the stability of Einstein-Maxwell theory near a Dirichlet 
wall. We focus on the case where the wall metric is Minkowski and the induced Maxwell 
potential vanishes on the wall. Far from the wall we take the spacetime to be asymptotically 
flat. For definiteness we work in 3+1 dimensions. 

We will use the term Dirichlet wall to denote any finite surface on which the induced 
metric is fixed. It remains an open problem whether the gravitational initial value problem 
remains well-defined at the non-linear level with such boundary conditions (though see 
[19-23]). We will not dwell on this issue here, but leave it for future work. 

In [4] we will show that Minkowski space with a flat Dirichlet wall is linearly stable to 
purely gravitational perturbations (though we will also find interesting instabilities in other 
contexts). A different sort of linearized analysis suggests that flat walls remain stable when 
exposed to simple non-perturbative processes involving neutral black holes. In analogy 
with electrostatics, positive mass sources near a Dirichlet wall must induce negative-mass 
image-sources behind the wall (see Fig. 1). Since the gravitational field of a negative mass 
is repulsive, any black hole is repelled from its image and thus from the wall itself. As 



a result, black holes in boxes should oscillate stably about the center of the box. This 
expectation will be verified explicitly in [4]. 

However, giving the the black hole an electric charge causes the image to acquire a 
charge of the opposite sign. The gravitational repulsion is then partly cancelled by the 
electrostatic attraction. In the extreme limit these forces cancel exactly cancel and further 
analysis is needed. 

We investigate this setting below, exploiting the fact that the cancelation of forces 
described above carries over to the full non-linear theory in the form of a family of static 
extremal black holes. As shown in section 2, this is precisely analogous to the cancellation 
of forces between extremal black holes. Indeed, our construction is the natural extension of 
that of Majumdar and Papapetrou [24, 25] to solutions with Dirichlet walls. The resulting 
static black holes provide a family of equal-energy solutions known as a moduli space. 
We then compute the moduli space metric i.e., the kinetic energy associated with motion 
through this space of solutions. Techniques for such computations were developed by Ferrel 
and Eardley to study colliding Reissner-Nordstrom back holes [26] and subsequently refined 
by Michelson and Strominger to study supersymmetric black holes in five-dimensional A f = 
1 supergravity [27]. The main subtlety in each method is the treatment near the horizon, 
where a naive calculation gives rise to divergences. Ferrell and Eardley motivate their 
regulator by considering lumps of extremal dust and noting that a singular limit recovers 
the Reissner-Nordstrom metric in Majumdar- Papapetrou coordinates (see (2.1) below). In 
contrast, Michelson and Strominger forgo this regulator and instead work directly with 
delta functions sources. Through clever manipulations they arrive at a finite expression for 
the moduli space metric without an explicit regulator. 

Applying analogous methods in our setting (see appendix A) yields a surprising result. 
Using h to parametrize the separation between the black hole and the Dirichlet wall, one 
finds the moduli space metric 


ds 2 


m 

~2 


dx 2 + dy 2 + 



9 

m 

i j 



(i.i) 


One sees that the coefficient of dh 2 becomes negative close to the wall, so that the black 
hole has negative kinetic energy. One might say that the center-of-mass motion of the black 
hole contains a ghost for small h. The computation is straightforward (but tedious) and is 
displayed in appendix A for those wishing to check the details. 

The result (1.1) is surprising enough that one should ask whether the techniques of 
[26] and [27] necessarily give the correct result in the present context. Indeed, though 
supersymmetry provides a strong constraint in the actual case studied in [27], we have 
found no proof in the literature fully justifying the treatment of divergences associated 
with the horizon. 

We therefore dedicate the majority of our work below to a complementary approach 
using Hamiltonian methods. Because they represent the total energy of a gravitational 
system as a boundary term at infinity, far from any possible divergence at the horizon, 
Hamiltonian methods are explicitly finite and gauge invariant. However, as we will see, 
they highlight the need for a boundary condition at the horizon to make the problem well- 
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defined. Below, we make the technically simple choice to cut off the near horizon region of 
the black hole throat by introducing second (topologically spherical) Dirichlet wall. 

Since the throat of an extremal black hole is infinitely deep, the original reasoning of 
Ferrell and Eardley suggests that moving the cutoff surface arbitrarily close to the horizon 
should not depend on the details of our boundary condition. But we will in fact find a moduli 
space metric that differs from (1.1). We therefore review this intuition more critically in 
section 2. In the end we will suggest - even in the absence of Dirichlet walls - that the 
moduli space approximation (also known as the adiabatic approximation) is ill-defined for 
classical extreme black holes due to the presence of long-lived excitations near the horizon; 
these are the excitations that were associated with turbulence in [28]. 

Interestingly, the additional term provided by our Dirichlet wall regulator fails to re¬ 
move the negative eigenvalue for black holes sufficiently close to the wall. Indeed, as we 
discuss below, the additional term contains a pole, and is associated with further pathologies 
of its own. Obtaining negative kinetic energies from both techniques suggests that this is 
the qualitatively correct physics and that kinetic energies can indeed become negative near 
the wall. This in turn indicates an instability to the black hole acquiring a large velocity. 
Unless some form of ghost condensation occurs, the Hamiltonian will be unbounded below. 

Such behavior would of course be forbidden in the absence of Dirichlet walls. For 
example, in flat space the positivity of kinetic energy on moduli space follows from the 
positive energy theorem of Gibbons and Hull [29]. By adapting the technology developed 
by Witten in [30], Gibbons and Hull show that E > \Q\ where E is the ADM energy and 
Q is the total electric charge. Since Q is also the energy of black holes at rest, the kinetic 
energy E — Q is non-negative. 

To investigate whether our Hamiltonian is bounded below, and also whether our neg¬ 
ative kinetic energies are in tension with the arguments of [29], we attempted to adapt 
the proof of [29, 30] to a spacetime with a Dirichlet wall. These attempts, however, were 
unsuccessful, leaving open the possibility that the theory is truly unstable. The difficulties 
with the proof may be connected to the fact that Dirichlet boundary conditions necessarily 
break supersymmetry [31] 1 . We briefly review the Witten style of proof and discuss the ob¬ 
struction in Appendix B. The essential point is that we are unable to construct the Green’s 
function necessary to solve = 0. 

The organization of the rest of the paper is as follows. In section 2 we construct a 
family of static extremal black holes parameterized by their mass and distance from the 
Dirichlet wall and derive a simple expression for their kinetic energy in the low velocity 
limit. Then, in section 3 we perform the integrals necessary to obtain the kinetic energy 
in closed from. Finally, in section 4 we discuss the physical significance of our results. We 
collect several technical results in further appendices. 


1 This fact is easily seen in the AdS context from the well-known result that holographic renormalization 
performed by imposing a finite-distance cut-off requires a cosmological-constant counter-term in the dual 
CFT (which would be forbidden if supersymmetry were preserved). We thank Will Donnelly and Joe 
Polchinski for discussions on this point. 
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2 Extremal black holes in the adiabatic limit 


This section develops the technology we will need to compute the kinetic energy of extreme 
black holes approaching a flat Dirichlet wall, where we use a second (roughly spherical) 
Dirichlet wall to regulate the calculation near the horizon and to impose a clean boundary 
condition. We first construct a moduli space of static solutions and then find an approximate 
solution for a slowly moving black hole. Hamiltonian methods allow the kinetic energy of 
this solution to be written in an integral form that will be explicitly evaluated later in 
section 3. 


2.1 The Moduli Space of Static Solutions 

Recall that the Einstein-Maxwell equations admit the Majumdar-Papapetrou (MP) family 
of solutions [24, 25] which take the form 

g a bdx a dx b = —ip~ 2 dt 2 + 0 2 (<fe 2 + dp 2 + p 2 d<p 2 ) A a dx a = 0~ l dt , (2.1) 


where 0 is any function satisfying 

B 2 ip = 0, (2.2) 

and d is the flat space (spatial) covariant derivative. The overbar denotes an exact solution 
at zero velocity and will be dropped for r / 0. 

We now impose boundary conditions on '0 at z = (1. Let E be a constant f-hypersurface 
in the half-space z > 0. The boundary <9E consists of two pieces, the surface z = 0 which 
we call <9Eo and the asymptotic boundary SE^. We require 

0 = 1 at <9Eo, and '0 = 1 + 0(r -1 ) near <9Eoo, (2-3) 

where r = y/p 2 + z 2 . These boundary conditions ensure that the solution satisfies 

• a Dirichlet boundary condition h ab = PAB, where Pab is the induced metric on the 
surface z = 0 (we use A, B for indices along the wall) and t)ab is the 2+1 Minkowski 
metric, 


• a conducting boundary condition (En = 0, B± = 0) on the vector held A a at 2 = 0 

• an asymptotically hat boundary condition far from the wall. 

If we allow -0 to be singular at a single point in E, we find the family of extremal black hole 
solutions 


0 


1 + 


m 


x — x+l 


m 

x — X - I ’ 


(2.4) 


where x± = ±hz. The integration constant m is the mass of the black hole (which is equal 
to the electric charge) and h parametrizes the separation between the black hole and the 
wall. The singular point (p = 0, z = h) is the location of the black hole horizon. 
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An interesting property of the solution (2.4) is that for h/m < 1/2 it contains trapped 
surfaces outside the black hole event horizon. This may be seen by defining the radial 
coordinate 

r+ := \x - x+\ = vV + ( z ~ h) 2 , (2.5) 

so that the area element y/ui on (topologically spherical) spacelike surfaces of constant re¬ 
takes the form 

Vw = r+0 = m+(l — r + + 0{r\) . (2.6) 

So for h/m < 1/2 the area of such spheres near the horizon decreases as we move outward. 
For a static solution this implies the existence of outer-trapped surfaces outside the black 
hole - a feature that cannot arise in asymptotically flat spacetimes satisfying the null energy 
condition without Dirichlet walls (see e.g. proposition 12.2.2 in [32]). This suggests that our 
Dirichlet wall may also lead to other features associated with violations of the null energy 
condition. We will confirm this suspicion below when we find negative kinetic energies. 
Interestingly, while we have found no trapped surfaces outside the black hole for h > m/2, 
the metric (1.1) implies we will nevertheless find that the kinetic energy becomes negative 
at h/m ~ 3.3, well away from h/m = 1/2, and the same will remain true with our Dirichlet 
regulator at r = a. 

2.2 Moving solutions and the great divide 

The geometric feature just described will turn out to force a division of our moduli space 
into two disconnected pieces when we add the additional Dirichlet boundary condition at 
r = a. But before discussing this “great divide” in detail, let us briefly recall just why such 
a regulator is needed at all. 

The point, of course, is that computing the kinetic energy requires approximate solu¬ 
tions in which the black hole has a small velocity v, so that the position of the black hole 
is time dependent. This entails promoting 0(a?) —» 0(f,x). But since since 0 diverges as 
x —> af_j_, simply replacing x± —> x± ±vt does not lead to a controlled expansion in v. Some 
regularization scheme is then needed to render 0 and its derivatives bounded so that we 
may meaningfully expand in v. 

We choose to introduce a radial cutoff a, calculate the energy for small v, and to then 
remove the cutoff by taking a —> 0. Let E be the t = 0 slice of our non-static spacetime and 
let dE be the boundary of E. The boundary dE now has three components (see Fig. 2): The 
asymptotic boundary dE^ far from the wall, the wall <9£o> and the radial cutoff surface 
dE a . We take dE a to be the codimension-two surface r+ = a and require the induced 
metric to be a (round) S 2 X M: 

ds 2 \ d z a = -dT 2 + R 2 dn 2 , (2.7) 

where dVt is the line element on the unit sphere and R is some fixed radius. The full set of 
boundary conditions on if) now reads 

= 1 at 8Eq , '0 = 1 + 0(r _1 ) near dEoo, and 0 = R/a at dE a . (2.8) 
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(a) 



d £ 0 


Figure 2. (a) A compactified representation of the Cauchy surface E and the three components 
of the boundary. dEo is the Dirichlet wall, 9E a is the cutoff surface (on which we also impose a 
Dirichlet boundary condition), and dE^ is spacelike infinity at which we require asymptotically 
flatness, (b) The black hole cutoff surface is a sphere centered at (z = h, p = 0) with radius a. 



Figure 3. A sketch of the polar coordinates (r±,0±) used in the text. 

Following [33], the solutions satisfying (2.8) can be approximated by the sequence of har¬ 
monic functions {ipk} where (see Fig. 3) 


R 


^ = i + (-- i )J2 ae+lB 


J +1 n (fc) ( Pe(cos{9 + )) _ ^ CQB(fl-)) 


J +1 


1=0 


(k) 

The multipole moments Bj, ; in the above expression are given by 

fc+l 

?(k) _ 




Here A = a/ (2 h) and the A p are given by 


(l + y/1 -4A 2 ) 1+P - (l - y/1 -4A 2 ) 
2!+P\/l - 4A 2 


Ap — 


i+p 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


We can now bound the error function A'Lfc := i)j ~~ V’fc- First note that vanishes 
on <9Eo,9£oo- By the maximum principle of the Laplace equation, |A?/>j,| must achieve its 
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maximum value on dT, a . In the limit of small a we have ~ A^ +1 and it is straightforward 
to derive the bound 

max (IA ipk I) < —A fc+1 + 0(X k+2 ) < aA (1_,5)(A;+1) (2.12) 

h 

for any k > 1. Here 5 may be arbitrarily small so long as a is appropriately large. Since 
A < 1/2, the sequence {ipk} converges uniformly to if:. 

In section 3 we will also need to approximate derivatives of ?/>. Arguments given in 
Appendix C show that for k > 1 

max (j<9q ... d in A^ k \j < P —. (2.13) 

Since multipole corrections carry no mass or charge, we can write the exact mass (and 
charge) m of the static solution as 

m= (R-a) lim B { Q k) . (2.14) 

k —^oo 


Inverting this expression and expanding for small a gives 

R = m + a(l- + 0(a 4 /h 4 ). 


(2.15) 


We may use (2.15) to eliminate R in a favor of m. and write ipi explicitly in the form 


4>i = 1 + {R — a) 


1 


1 + A + 


A 2 


1 — A 2 


1 


1 


-) — a ( A + 

r + r_ 


A 3 


(1-A 2 ) 2 


1 


= 1 + m (-) — 

r_|_ r_ 


ma 3 ( cos (0+) cos(0_) 

J2hy ' " + 


+ . . . . 


cos(0 + ) cos (6-) 


(2.16) 


where x± = ±hz, (r-t, 6±) are polar coordinates centered at x± as in hgure 3, and the ... 
in the second line represent terms which are subleading as a —> 0. 

However, (2.15) turns out to have a more fundamental implication. For h > m/2 with 
a small and positive, it is clear that R > m as one expects. But for h < m/2 one finds 
R < m. This is another manifestation of the geometric result discussed around (2.6). But 
since R is fixed by our choice of induced metric at r = a - i.e., as a boundary condition - 
this requires us to use different Dirichlet regulators on each side of h = m/2. We are unable 
to regulate the entire moduli space at once. This is the “great divide” mentioned earlier. It 
suggests that the two pieces may fail to join smoothly at h = m/2 even when the regulator 
is removed. 

It is natural to ask if one could do better with another boundary condition. Rather 
than investigate specific alternatives, we conclude this section by estimating the dependence 
of our final results on the particular choice of boundary condition at r = a. One might 
expect that the infinitely deep throat of the extremal black hole would remove sensitivity 
to the boundary condition as a —> 0, but this intuition can be tested as follows. Let At be 
the minimum time, as measured at <9Eo, for a light signal to make a round trip from <9£o 
to dTi a and back. This is the timescale on which we expect the boundary condition at <9£ a 
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to become relevant to the motion of the black hole. We can then calculate Ah, the change 
in h over a time At. The result is 

Ah = vAt = v f = v f drip 2 = mv (— + 21og(/i/a) N ) + ..., (2-17) 

Jh r J h V a ) 

where ... denote terms that do not diverge as a — > 0. 

In order to obtain a controlled expansion in which the metric backreaction is small, it 
will be important below to take v —> 0 such that ip ~ vdip ~ mv/a 2 does not become large; 
say that mv/a 2 remains less than some constant c. Then (2.17) gives 


Ah 


ma 


< 


c. 


(2.18) 


This Ah is the distance over which causality protects finite motions on moduli space from 
dependence on the details of the boundary condition at r = a. And (2.18) shows that 
Ah —> 0 as a —> 0. So in this limit causality turns out to provide no protection at all. 

This suggests that the adiabatic approximation is not in fact well-defined for classical 
extreme black holes. Instead, it will depend on a choice of regulator, perhaps implemented 
as a boundary condition at the horizon. 

Nevertheless, we would like to investigate some specific case in more detail in the hopes 
of extracting useful qualitative features. To this end, we will carry a detailed calculation 
with the particular Dirichlet regulator at r = a described above. We will compute 


T = lim 

a —>-0 



E — m\ 

) ’ 


(2.19) 


were E is the total energy and m is the energy of the static solution. In addition to an 
ADM-like boundary term at infinity, we expect E to have contribution from the Brown- 
York stress tensor [3] at the wall. At least with our choice of regulator, v 2 T is the leading 
order approximation to the kinetic energy of the black hole. 


2.3 The adiabatic expansion 

In order to study the kinetic energy v 2 T, we construct approximate solutions for small 
v. We wish to work in the adiabatic limit, in which we can associate our solution with a 
trajectory through the moduli space of (regulated) static solutions. In practice this means 
that we replace the boundary condition (2.8) with 

ip = 1 at BE o, ijj = 1 + 0(r -1 ) near SSqo, and ip = R/a(t ) at dE a (t). (2.20) 


We take the boundary dT, a ^ to be the surface r + (f) = \Jp 2 + [z — h(t)) with a(t ) is defined 
in terms of h(t) by using (2.14) and requiring that R , m be time independent. At time t = 0 
we let h = —v, where an over-dot denotes a t-derivative. We find 


mav 

h{2h — m) 


+ 0(a 4 ). 


( 2 . 21 ) 


At each time t we take ip to solve the Laplace equation (2.2) with boundary conditions (2.20). 
This defines ip(x,t). 





We will need the time derivative ip to compute the kinetic energy below. Since ip 
is a solution to the Laplace equation, it is uniquely determined by its boundary data. 
Using (2.20), we can write down boundary data for ip: 

ip = 0 at <9S 0 , ip = 0(r -1 )near dt,^, ip = ad a - f+d r+ ip\ r+=a at <9£ a ( t) . 

(2.22) 

The first term at dT , a m is the time derivative of the boundary condition on ip , while the 
second term accounts for the fact that dT, a is moving through space. The partial derivative 
d r+ is taken after expressing ip in terms of the coordinates (r_|_,0_|_) from figure 3 using the 
relations 

r_ = y/(r+ sin(0+)) 2 + (r+ cos(6>+) + 2 h) 2 , tan(6L) = — r+ 9 , • (2.23) 

r + cos (0+) + Ih 

We may then find r + by solving 

( r + pt ) 2 sin 2 (0 + ) + [r+(t) cos (6+) + vt ] 2 - apt ) 2 ) = 0, (2.24) 

and setting r + = a. This yields 

r+\ r+ =a = a-vcos{6 + ). (2.25) 

Below, we will approximate ip by a function that is precisely harmonic but only approxi¬ 
mately satisfies the boundary conditions (2.22). We will construct such a function at the 
end of section 2.4. 

For i) / 0 the full spacetime metric and Maxwell potential contain additional compo¬ 
nents not found in (2.1). A general ansatz is 

ds 2 = —{ip + 5ip + 5ipt)~ 2 dt 2 + 2 ip~ 2 R ■ dxdt + {ip + 5ip) 2 {5ij + <7ij)dx l dx\ 

A a dx a = — {ip + dip + 5ipA_)~ 1 dt + {P — ip~ l R) ■ dx, (2.26) 


where a t j is traceless. We also require the solution to be invariant under the combined 
transformation t —>• —t and v —> —v, and to reduce to (2.1) when v = 0. This means that 

R,P~v, 5ip,5ipt,SipA,crij ~ v 2 . (2.27) 


The form of the ansatz was chosen so that R transforms simply under the coordinate 
transformation t —> t + n{x), namely as R —> R + dfi. 

By calculating the induced metric and gauge held on the codimension-one surface 
r + (t) = apt ) we see that our Dirichlet boundary conditions require 

R = ^vip 4 {f ■ z)r, P =-vip 4 {r ■ z)r (2.28) 


at this surface. The situation is simpler at z = 0 where we need only require R, P to be 
normal to the wall. Asymptotic flatness further requires that the metric and gauge held 
satisfy 


1 

Qab Vo,b ^ — 



(2.29) 
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Boundary conditions on the second order quantities 5'il>,5'ipt,3' l pA,c r ij at <9Xo and <9E a can 
be also worked out from our Dirichlet boundary conditions, though they will not be needed 
below. 


2.4 A simple expression for the kinetic energy 

We may now compute the kinetic energy v 2 T. As noted in 2.2, in addition to the familiar 
ADM-like boundary term at infinity, the energy term in T should receive a contribution 
from the Brown-York stress tensor at the wall. One could proceed by computing both 
contributions and manipulating the resulting expressions. However, we will instead use the 
Landau-Lifshitz technique [34] to express the total energy as a bulk integral of the graviton 
stress tensor (see also the closely related approaches of [35, 36]). 

The basic idea of the Landau-Lifshtiz technique is to note that, about any background 
solution, the linearized Einstein equations define a tensor that is an identically- 

conserved tensor with respect to the background covariant derivative V a ; i.e. V a i7^ ab = 0 
for any linearized fields (whether or not they satisfy the linearized equation of motion). So 
if the background has a timelike Killing field £“, integrating £ a n b H^ over a Cauchy surface 
with normal n b defines a conserved quantity even when the fields are not small. In other 
words, we may replace what were formally tangent vectors to the space of solutions with 
the finite differences between the actual fields in any solution and fields in the background. 
One then notes that £ a n b H^ is a total derivative, so as expected this energy can be written 
as a boundary term 2 . A major advantage of this approach is that it does not require us to 
solve the second order equations of motion. In fact, as in [26, 27], we will see that it is not 
even necessary to solve the linearized equations of motion. 

The same idea can be applied to the electric charge. So we may compute the kinetic 
energy (2.19) using the conserved current 


Ja ■ = 


r(l) 


b 77(1) 


CH J - -V°F f 




ab 


In addition, £ is the background timelike Killing vector and 


H ab := G ab - 8 


(2.30) 


(2.31) 


where k = 87T, and £ = —dt is the background (v = 0) timelike Killing vector. The 
superscript (1) indicates that we keep only terms linear in the metric and Maxwell fields 
(though these terms can be non-linear in the velocity v). A straightforward calculation 
gives 


kJ = -dtZ, kJ q = ■ Z, 

Zi = 2^di(^ 2 5^ A ) + T^djCTij + 2- 2(2.32) 

2 At least to the order used here, the equivalence of this approach with other standard definitions fol¬ 
lows from e.g. [37] which, as shown explicitly in [38], implies that when perturbing around a stationary 
background, the lowest-order change in the Hamiltonian is the sum of the above Landau-Lifshitz term and 
the change an appropriate horizon area. The latter term vanishes in our context since the region between 
infinity and our wall at r = 0. is horizon-free. It then follows that for us the lowest-order change occurs at 
what we have called second order, and that this is properly computed by Landau-Lifhshitz techniques. 
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where we have used the background equations of motion 3 2 ip = 0 and the linearized equa¬ 
tions of motion. One may then verify directly that 


V“ J a = 0, 


(2.33) 


where V is the background covariant derivative associated with setting v = 0. 

The associated charge is 

Q ■= ~ j V^u a J a , (2.34) 

where u a and 7 are the background unit normal and induced metric on the X (which was 
defined above to be the surface t = 0). Since (2.19) may be written 


T = lim 

a —>-0 


O 

lim 

v —>0 V A 


(2.35) 


we need only compute Q to second order in v. It is straightforward to check explicitly that 
Q is conserved to this order using 

Q = d t (^-[ d 2 y n- z\ = — f d 2 yn-J, (2.36) 

V® 71- Jot. J Jan 


where n is the unit normal to <9X. By time reversal symmetry J can only contain odd 
powers of v, but the terms of order v vanish by the linearized equations of motion. 

We now derive a convenient expression for computing Q. We begin with the equality 

H 00 + -V 6 % = 0 , (2.37) 

V 

which follows from the Einstein-Maxwell equations. Taking (2.37) and putting all of the 
linear terms on the left hand side gives 


Z 


3 ip 2 \d x P \ 2 4(3 x P) ■ 0 x R) 9|3 x R\ 2 

y) 2 + ip 6 ip 7 + 4^ 8 


di 


+ 


2 ip lr tpRi + d[i(ip A Rj]Rj) 


ip 4 


(2.38) 


To simplify this expression, we must gain greater control over P, R. Using the ansatz (2.26) 
and expanding the equations of motion to linear order yields [26] 


d x P d X R 

ip 2 ip 3 

d x P 3d x R 

ip 3 4 ip 4 


S X K + 34-1 

34 » 2 , 


where 4>i, 4*2 are integration functions and K is defined by 


d ■ K = —ip B 2 K = 0. 


(2.39) 

(2.40) 
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These equations are equivalent to the gravitational momentum constraint 13 = 8irn a Tfi M 
and Maxwell equation X7 a F aj = 0, which are the only non-trivial equations of motion at 
linear order in v. 

Noting that each of K, R, P are invariant under (p —> <fi + const and <f> — q i, their 
curls dxK,dxR,dxP must be invariant under cp —> cp + const but odd under <p —>• —(f>. 
As a result, only the (^-components of these curls can be non-zero. But the symmetry also 
requires dq, dq to be independent of cp, so the ^-components of the gradients Sdq, <9d>2 must 
vanish. 

It follows that the curl terms and gradient terms in (2.39) vanish separately. As a 
result, we must have 


d x R = -4ip 3 d x K 

B x P = -3tp 2 d x K. (2.41) 

Combining (2.38) and (2.41) (along with the definition ip = —d ■ K) we obtain 


d-Z = —3ip 2 


d-K + 


dxK 


+ di 


2 ^Rj 




(2.42) 


Inserting this expression into (2.34) gives 


Q = — / d 3 xj? 


d-K) + 


d x K 


+ k Jqy; " ' ^ ip 


ni ( + d[i(i/> 4 Rj]Rj) ) , (2.43) 


where n is the outward pointing flat space, unit normal to dT,. The second term in the 
surface integral vanishes by the Dirichlet boundary condition which requires that R is 
normal to the <9£. The first term vanishes everywhere except on the cutoff surface (for 
which ip ^ 0). We may thus write 


Q = -— / d s xr 

8vr Jr 


d-K + 


dxK 


+ 


3 vR? 
87 r a 


ldY.a 


dO, xp cos(0 + ), 


(2.44) 


without having to solve the linearized constraints. 

Computing Q analytically requires us to write down closed form approximations of 
d ■ K and | d X K\ as we did for vp above. Such approximations are obtained by replacing 
K with 


K\ 





/ 1 1 \ mva 3 / cos(0 + ) 

\r+ r_ J 4 h 2 \ r\ 


cos(6L)\ 

) 


z. 


(2.45) 


An explicit calculation reveals that ipi := ~d • K\ satisfies the boundary conditions (2.20) 
up to terms of order 0(a). In appendix D, we also show that \d x (K — K\)\ ~ 0(a). As 
a result, we will see in section 3 below that replacing (ip,K) —> (ipi,Ki) gives us sufficient 
accuracy to compute the kinetic energy (2.19) in the a —> 0 limit. 
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3 Computing the kinetic energy 


As a warmup we now calculate the kinetic energy without the wall at z = 0 and obtain the 
expected answer mv 2 /2. In the absence of the wall the functions ip,ip,K are precisely 


, m ■ mvcos(d) -> mr. 

ip = l + —, ip = - K= -z 

r f 


(3.1) 


where r = p 2 + (z — vt) 2 and m = R — a so that ipdT, a = R/a. Inserting these expressions 
into (2.44) gives 


Q = -- 


dr 


m 2 v 2 (m + r) 2 3 v(m + a) 


+ 


d(cos(8)) 


mv cos 2 {0) 


mv 2 (R 3 \ mv 2 (R 3 \ mv 2 
—>— 1 I H - x - ( —5- 1 = 


2 ’ 


(3.2) 


as expected. Note that both terms in (2.44) diverge as a —> 0, though the divergences 
cancel and leave a finite answer. Interestingly, for this case our Dirichlet regulator at r = a 
gives precisely the same result as the techniques of [26] and [27]. 

We now return to the case with the wall. First we consider the boundary term in (2.44) 
with K replaced by K\. From (2.45) we find ipi := -d-K\. Performing the angular integral 
then gives 


^- [ dnipcos(0+) 

OTTCL J r +— a 


mv 2 (R 3 m?(m + h) \ 

2 \ a 3 4/i 3 (/r — m/2) J 


+ O(a). 


(3.3) 


As a simple check on this result, note that if we take h —> oo we recover the boundary 
term that appeared in (3.2). As with ipk we could define a series of functions which 
include higher multipole corrections. Power counting shows that replacing ip\ with ^ with 
k > 2 would not change any of the terms shown above. So the approximation K\ suffices 
to compute this term in the limit a —> 0. 

Evaluating the bulk term is somewhat more involved. After again replacing K by 
K i, we make the decompositions Vh = V’o + an d K\ = Kq + 5K. Let us first study 
terms involving 5ip, 5K. All but two of these terms vanish by simple power counting. The 
remaining terms are 


3 

87T 

3 

87T 


d 6 x' 0o 2 (d ■ A'o)(<9 • SK) + (d x A'o) • (d x 5K) 


d 3 x2il’o dip 


(d-K 0 ) 2 + \dxK 0 \ 2 = 0(a), 


m 4 (m + 4 h)v 2 
8h 3 (h — m/2) 


+ 0(a) 


(3.4) 


where we have now evaluated these terms by using power counting to show that, in the 
first line, only the leading order term contributes as a —> 0 and, in the second line, only the 
spherically-symmetric part of the term in square brackets can contribute. The second line 
then vanishes as a —> 0 since d'lp contains no monopole term. For the same reason, further 
improving the approximation by passing to K n for n > 1 gives no effect as a —> 0. 
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It remains to calculate the terms involving only ip o and Kq. Here it is convenient to 
make contact with the technique developed in [26, 27]. Let us first rewrite tpQ. Kq as 


^o = l+ V —, 


A=± 


ta 


*0= E 

A=± 


m A 

r A 


-v A , 


(3.5) 


where rn± = ±m and v± = =p vz. We then introduce the derivative operators d± which act 
on r± = \x — x±\ by taking gradients with respect to x±. We will treat x± as independent 
parameters until the end of the calculation at which time we set x± = ±hz. Using this new 
notation we can write 


(V • Kq ) 2 + |V X Kq 


■ v b 9a ip ■ d B ip - VA ■ d B ipv B ■ dAip + VA • 9a^v b ■ d B ip) 

A.B 

v 2 (d + ip ■ d+ip + d-ip ■ d-'ip — 2 d + ip ■ d-ip). (3.6) 


After some straightforward algebra we find that this contribution to Q can be written 



where 


V = d + ■ d + + d- ■ d- — 2d + ■ d- . 


(3.8) 


Here we have exploited the fact that T>ip^ can be written as a total divergence. This allowed 
us to separate the integral on the first line of (3.7) into a boundary term on dT ,o U clSoo 
that remains bounded as a —> 0 and a boundary term at dT, a that diverges. The boundary 
term on evaluates to the first term on the second line. The boundary term <9£o U dT,^ 
can then be written as a bulk integral over the entire region z > 0 (including the region 
r < a) treating the integrand T>ipQ as an appropriate distribution at r = 0. This term is 
evaluated in Appendix A, which also explains the connection of this term to the effective 
action of [26] or [27]. 3 The result is 

- aa [ Vipl = T ^~ (1 - U(m/h )), with U{x) = \(6x + 2x 2 + x 3 ). (3.9) 

'j^ 7r Jz >0 * * 


Combining (3.3), (3.4), and (3.7), the divergent i? 3 /a 3 terms cancel to yield 
1 If 6 m 2 m 2 m 3 \ 3m 3 


Q = mv A 


6 m 2m" m" 

h h 2 h 3 ) 8h 2 (h — m/2) 


+ 0(a). 


(3.10) 


The hnal term inside the brackets (the term with the pole at h = m/2) is the sole discrep¬ 
ancy between our Dirichlet-regulated calculation and the result one would obtain using the 
techniques of [26, 27]. Working backward, one may show that this term originates from the 
dipole correction to ip appearing in (2.16). This fact strongly suggests that our hnal result 
for the kinetic energy does in fact depend on our choice of regulator. 

3 cf. Eq. (2.10) in [27] after setting A = 0. Note that in [27] the authors work in five spacetime dimensions, 
whereas we are working in four. See Appendix A for additional details. 
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4 Discussion 


Our work above computed the metric ou the moduli space for extreme black holes in the 
presence of a flat Einstein-Maxwell Dirichlet wall using two distinct methods of dealing 
with divergences at the horizon. The results indicate two important lessons. The first 
stems from the negative eigenvalue of the metric for black holes close to the wall. Since it 
appears using either of our methods, the associated negative kinetic energies appear to be 
a robust feature of the physics, with corresponding implications for stability. 

Indeed, let us imagine that we place an extremal black hole in the region with a negative 
eigenvalue. Then the negative kinetic energy means that conservation of energy allows any 
process emitting positive-energy gravitational radiation so long as the black hole acquires 
a non-zero velocity v. Furthermore, the process can continue (with corresponding increases 
in velocity) at least until v becomes so large that the moduli space approximation breaks 
down. Referring to this component of velocity as a ghost, the system must either exhibit 
some form of ghost condensation (which would not be visible in perturbation theory) or the 
Hamiltonian will be unbounded below. We speculate that the latter may in fact be true. 
In this direction, we note in appendix B that standard techniques for proving a positive 
energy theorem [30] break down when we include a Dirichlet wall. The same comments 
apply to attempts to use techniques from [29] to prove E > Q. 

The second lesson stems from the fact that our two method gave results that differ in 
detail. Indeed, adding an explicit cut-off at r = a led to a new term in the moduli space 
metric beyond the result (1.1) that would be obtained using the techniques of [26] or [27]. 
With a Dirichlet boundary condition at this cut-off surface, we found a particular such term 
containing a pole at h = m/2 that remains after taking a —> 0. We also saw that the origin 
of this term was the dipole correction which must be made to to satisfy the Dirichlet 
condition at r = a for finite a. Since dipole corrections should be expected to arise from any 
other choice of boundary condition at r = a, this suggests that other boundary conditions 
again lead to a —> 0 moduli space metrics that differ from (1.1), and that the result may 
depend in detail on the boundary condition chosen. Our moduli space of extreme black 
holes would then be ill-defined. 

It is natural to ask whether this regulator-dependence arises only in the presence of 
the Dirichlet wall <9Eo. We have therefore also used our Dirichlet regulator to calculate the 
moduli space metric for the two black hole scattering problem originally treated in [26]. 
For that case we find the Dirichlet-regulated metric to again contain an additional term 
beyond those appearing in [26]. The new term is precisely the one appearing in (3.10) with 
m replaced by —m. The dependence of the moduli space metric on regulating boundary 
conditions thus appears to be a general property of moduli spaces for classical extreme 
black holes. 

Indeed, we see no a priori reason to expect a good adiabatic approximation even as 
v —> 0. The point is that the lifetime of the longest-lived quasinormal mode of a classical 
black hole turns out to diverge as the black hole becomes extreme. It is precisely this 
phenomenon that was associated with turbulence in [28]. And for any finite v process, 
it suggests that a sufficiently extreme black hole can be significantly perturbed from any 
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nominal “ground state.” It would thus be very interesting to better determine if (and 
when) the moduli space calculations of [26] correctly approximate the kinetic energy of an 
interacting pair of slowly-moving extreme black holes. While numerical studies (see e.g. 
[39]) have not yet reached the level where they can be usefully compared with such results, 
we may hope that they will do so in the near future. 

It remains to address the success of [27] in comparing its classical moduli space with 
string-theoretic calculations. If the classical moduli space is ill-defined, why should there 
be any such agreement? Recall that the goal of [27] work was to compare a classical black 
hole calculation with the moduli space for a supersymmetric quantum system. Quantum 
mechanics should render the black hole spectrum discrete, so that an adiabatic approxi¬ 
mation will in fact hold for sufficiently slow processes 4 . Thus the desired quantum system 
should indeed admit a well-defined moduli space. One may thus view the success of [27] 
as suggesting that requiring the moduli space to maintain supersymmetry is sufficient to 
extract the associated metric from a classical calculation. 

Having argued that the a = 0 moduli space is ill-defined, we pause to note that space- 
times with finite cutoff parameter a are interesting solutions which form a moduli space in 
their own right (for a given sign of h — m/2). Taking a to be small but finite we may inter¬ 
pret (3.10) as the kinetic energy of a solution with two Dirichlet surfaces, one flat and one 
spherical. Our results clearly show that such boundary conditions admit negative kinetic 
energy solutions with E < Q. 
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A The Bulk Integral 

This appendix derives (1.1), which gives the moduli space metric for extreme Einstein- 
Maxwell black holes near a Dirichlet wall as computed using techniques of [26, 27]. This 
method focuses on the effective action 

5eff= I(b / d4x ^ R - FabFab )+ 

J ( U a - A a )dx a j}j . (A.l) 

4 Since the spectrum will generally have exponentially small spacing, this moduli space metric will govern 
only exponentially slow motions. The notable exception occurs in supersymmetric contexts where the gap 
above the highly-degenerate ground state is only polynomially-small. 


5 , 


source — , 

47T 
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The term 5 SOU rce couples extremal point particles of of mass and charge m with four velocity 
u to the fields. In [26], the source terms were derived by writing down the action for an 
extremal dust and taking the dust distribution to be a delta function. In [27] these terms 
are reinterpreted as compensating for the fact that the field equations are not satisfied at 
the timelike singularity behind the black hole horizon. 

We are to compute this effective action for the ansatz 

ds 2 = —ipQ 2 dt 2 + 2 i/i~ 2 R ■ dxdt + V>o(<% + aij)dx l dx ^, 

A a dx a = —^Q l dt + (P — V ; o" 1 -R) ' dx , (A.2) 


with i/jq as given as in (3.5). After a long but straightforward calculation one finds 

Seff = - 3^- £ f dt(6 l] 5 kl + SiSi - 5\5{)v k A v l B d Ai d BjX 

n A =±,B=± J 


=- T>\ , 

32vr A ’ 


(A.3) 


where we have used the linearized constraints to eliminate R, P from the action, 


X = 


[ 

J 2>0 


(A.4) 


and 


V = d + ■ d + + d- ■ <9_ — 2d + ■ <9_ . 


(A.5) 


This is a four-dimensional analogue of the 5-dimensional result of [27]. 

The definition (A.5) agrees with that of (3.8) from the main text. While y diverges, 
T>x may be evaluated by first acting with V on the integrand ipQ, treating the result as a 
distribution, and finally performing the integral over all 2 > 0. With this understanding 
we see that T>x agrees with the term evaluated in (3.9). 

To begin the computation of Dy, let 


where 


X = do + mil + rn 2 h + m 3 h + m 4 / 4 , 


Iq= ux 
J Z>0 


h = 4 / d 6 x 

J 2>0 


1 1 

r + r_ 


I 2 = 6 / d 6 x 
J z>0 


1 1 2 

2 T 2 

r± r_ r + r_ 


I 3 = 4 / d A x 
J z> 0 


1 1 3 

’>■ *7" + 


3 2 1 2 2 

y* ' J ry* U ry> ^ ry> ^ ry* ^ ry» 


I4 = / ax 

J z> 0 


1 | 1 | 6 4 4 

A A 9 O’ Q - Q 

r + r+ri r+ri 


(A.6) 

(A.7) 

(A.8) 

(A.9) 

(A.10) 

(A.ll) 
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As advertised above, we evaluate VI n by acting with P on the integrands (A. 8), (A. 9), 
(A. 10) and (A. 11) and calculating the integrals afterwards. We take the the black hole 
and its image to he on the z-axis in positions which we denote by 2 ;+ > 0 and < 0, 
respectively, and after taking the desired derivatives we set = — z+ as required by our 
boundary conditions. 

In the spirit of [26, 27] we integrate over the entire region 2 ; > 0. as a result, we will 
often use the identity dlrZ 1 = —47 t5^ 3 \x — x+). We now collect the useful results 


r + = r_, 


d z 


nz + 

31+2’ 

r + 


d z r 


—n 


nz + 

31 + 2 ’ 

r + 


at z = 0, 


(A.12) 


d 2 xr+ n = 


2vr 


zl n for n > 2, 


/*= 0 fa - 2 ) 

d 2 rZ 2 = 2rZ 4 + 2rZ 1 h^ 3 ^(x — xJ) 


,q _2 —4 7T 

drxrrr_ = 


' 2>0 


2^3 ' 


(A.13) 
(A. 14) 
(A.15) 


We emphasize that identities (A.12) and (A.15) hold under the assumptions that x+ = 
— X- = z+z, which we impose after computing all the derivatives. We also mention that 
identities (A. 14), (A.15) are useful to compute the term V f r^ 2 rZ 2 . Finally, we shall make 
extensive use of the fact that x_, the position of the image, lies outside of the domain of 
integration of the integrals (A. 8), (A. 9), (A.10) and (A.11), so volume integrals involving 
<5( 3 )(if — X-) vanish. 


A.l Terms with m 2 

We calculate now the quantity P/ 2 , where 

I 2 = 6 / d 3 x 
J z> 0 



2 

r+r- 


■ term r, 2 


dl 





d+ ■ d- f r + 2 = 0 
J z> 0 


2vr 

-+ 


Hence, 


■ term r 



2vr 

z+' 


d 2 _ 




o 2 -2 

o r_ 



d z r 


-2 



2ir 

z+ 


(A.16) 

(A. 17) 
(A.18) 

(A.19) 

(A.20) 
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(A.21) 


d 2 + I rZ 2 = 0 
J z> 0 

d + ■ d- f rZ 2 = 0 
J z> 0 


Thus, 


V 


r~ 2 = 


<z> 0 


term r, V 1 


d. 


+ '+ 


r i _ 1 r _ 1 = — 7 


2vr 

«+’ 

47T 


2vr 


'z>0 

where we have used in the last step that |T + — x_| = z+ — z~- = 2z+. 


d: 


r , V 1 = 


' z>0 


!z> 0 


r~ 1 d 2 rZ 1 = o 


since x_ lies outside of the domain of integration. 


(A.22) 
(A.23) 


(A.24) 


(A.25) 


d + -d- 


' z>0 




Vr , 1 • Vr 


-l 


' z>0 



= 2 + 



So 


and 



47T 

5 


X>/ 2 


487T 

2 + ’ 


A.2 Terms with m 3 

Now we calculate P/ 3 , where 


h 


= 4 


[ d 3 x 
Jz> 0 

' 1 

1 3 

I 

3 

L r + 

T 1 9 

ry* ^ tyt 

2 

r+r- 


— (A.26) 
^+ 


(A.27) 
(A.28) 


■ term r , 3 


Hence 


■ term r 




d+ ■ d- f ?’ + 3 = 0 

J z>0 


V 






(A.29) 
(A.30) 
(A.31) 

(A.32) 

(A.33) 
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(A.34) 


/ r_ 3 = 0 


<9+ ■ 8 - / ?’_ 3 = 0 


(A.35) 


Hence 


V / r~ 6 = 


(A.36) 


term r, V 2 


Bi / rTV - 2 = - 


|f+ - ,f_| 2 z 2 


(A.37) 


£ / rTV- 2 


A 1 9 2 r " 2 


/ ?+<9 2 r_ — r_^<9 2 r_ + / r_d 2 r J[ 

' z=0 Jz >0 

47T f _= 


|x+ — x_| 


<9 + • <9_ / r,,, 1 ?-. 2 = / V?’ + 1 Vr_ 2 

J z>0 J z> 0 


r_ 2 d z r_L_ l — / r_ 2 d 2 r , 1 


rw-^12 “ z + / r 4 

F+-®-l 4 2 =o 


(A.38) 


(A.39) 


So we have 


£> / r -l r - 2 = - — 

,/^n + " 3z 2 


(A.40) 


term r , :2 r 1 


32 / _-2 -1 


Q 2 r A 2 r_ l = — / r_ 1 0sr + 2 — r + 2 9 2 r_ 1 = —3z + / r_,_ 5 = —tt 


5 _ 

r + “ .2 


(A.41) 


^ / r; 2 rl 1 = 0 


(A.42) 


9 + • <9_ / r , 2 r_ = / Vr, • Vr_ = — / r + 2 d z r_ =z + r, b = —tt- (A.43) 
Jz>0 Jz> o Jz=o Jz =o 3z, 
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Thus, 

and we have 

A.3 Terms with m 4 

Finally, we compute VI 4 , where 

r 

A 


V 


>z> 0 


r -2 r -l _ _ 

r + r ~ ~ Q ~2 ’ 


107T 


w 3 = 


167T 

4T 


(A.44) 
(A.45) 


h = I dr x 
J z> 0 


116 4 4 

A “I - 4 ~0 o~ Q - 

r + r± r+r_ r+rl 




terms r , 4 and r 4 


The computation of these terms is closely analogous to the one of T>r .|_ 3 . We obtain 


v I T = -% 

lz> 0 Z+ 


Vj 

Jz>o 24 - 


term r, 4 r 3 


3 


47T 


7T 


d+ Jz> 0 T+r ~ |f+-f-| 3 24 


(A.46) 
(A.47) 

(A.48) 


^ / r + 1 r_ 3 

'z> 0 


4 


f, 1 d 2 r_ 3 


J z> 0 

47T 

|.T_|_ — X— | 3 
47T 

|.T + — ,T _| 3 
37T 

24 


/ (r + 1 d 2 r_ 3 - r_ 3 d 2 r + 4 ) 

J z=0 


1 z =0 

+ 4z + / r. 
Jz=0 


-6 


(A.49) 


Thus, 


<9 + • 4 / r + 4 r_ 3 = / Vr + 4 Vr_ 3 

«/ £>0 J z>0 

47T f 


\x+-x- 1 3 4=0 

47T I 


r_ 3 <9 z r, 1 


\x+ — X- 


= 0 


V I r , 1 r 3 = - 5 -. 


\3~ Z + r 
Jz =0 


7T 


-6 


<z> 0 


+ _ ~3 


(A.50) 
(A.51) 
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term r, 3 r 1 


Si 


' z> 0 


r-V 1 = - 


f {r_ l 8 z r + 3 - r + 3 d z r_ l ) = -4 z + f r + 6 = - ^ 
lz =o Jz =o 2+ 


Si I r+Vz 1 = 0 


>z> 0 


7r 


<9+-3- / r + d r_ =- / r + 6 d z r_ = z + J r + = ^ 3 "- 


' 2 >o 


' 2=0 


Thus, 


P / r+Vz 1 = -Tr 

J 2>0 


37T 


term r , 2 r 2 


(A.52) 
(A.53) 
(A.54) 

(A.55) 




-2 -2 

r, r_ = — 


' 2>0 


I (r_ d z r + - r + 8 z r_ )+ r + 8 r~_ 

Z =0 JZ >0 


= -4 

7r 

= -la" 


[ r + 6 + 2 / rZ 2 rI 4 

J z =0 «/ z>0 


(A.56) 


Here we have used identities (A.14) and (A.15). These identities may also be used to show 

(A.57) 


o2/ -2-2 / —2 o2 —2 r, / -2 -4 7r 

8-1 r + r_ = j r + 8 r_ =2 / r + r_ = 


' 2>0 


' 2>0 


' 2>0 


and 


<9+ • <9_ f r + 2 r_ 2 = — f r + 2 d z r_ 2 — f r + 2 d 2 r_' 
J z> 0 J z =0 J z> 0 


= 2 z 4 


= 0 . 


/ r + 6 — 2 / r + 2 r_ 4 

' 2=0 Jz >0 


Hence, 


P 


r^ 2 rl 2 = 0, 


' 2>0 


§77 

= Z3~ ‘ 


and 


A.4 Final result 

The action of the operator P on the various integrals J n was found to be 

P / 0 = 0 P / 1 = -167T P/ 2 = 4 f P / 3 = ^J Vh = % 

where we have used z_|_ = h. Inserting these results into the expression 

v 2 v 2 

- T>X = - v Y" m n I n , 

327r 39 tt n ’ 


32tt 


(A.58) 

(A.59) 
(A.60) 


(A.61) 

(A.62) 


gives ( 1 . 1 ). 
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B Difficulties of proving a positive energy theorem 


We now review Witten’s method of proving a positive energy theorem and show why this 
technique appears not to extend to settings with Dirichlet walls. 

Witten’s argument [30] consists of two parts. One first supposes that the spacetime 
admits a Cauchy surface £ and a spacetime spinor satisfying 


^6 = 0 , e = eoo + 0(r x ) on £, 


(B.l) 


where Coo is a constant spinor, V a are the spacetime covariant derivatives, 7 “ are the Clifford 
algebra generators satisfying 'y a 'y b = — 2 g ab , and ^ = 7 with an index i denoting 

a corresponding object projected into the surface £. A calculation then shows 


(eooeL)^ > _/jVe| 2 > 0, 


(B.2) 


where E is the energy associated with a future-directed asymptotic Killing held proportional 
to CooT^oc- The inequality (B.2) follows after rearranging 


[ e t y 2 e = 0, 

is 


(B.3) 


and assuming the dominant energy condition. 

It is then argued that a solution to (B.l) always exists for asymptotically hat space- 
times. This is shown by decomposing e = e\ + € 2 , where is a trial function of the form 
£i = £oo + 0(r~ l ), so that 62 ~ 0(r _1 ) for large r. Then 62 can be written formally as 


e 2 (x) = y dyG{x,y)[^7ei}(y), 
were G is a Green’s function G satisfying 

rfG{x,y) = 5(x, y). 


(B.4) 


(B.5) 


In other words, G = D , where D is the differential operator ^ restricted to spinors that 
vanish as r —> 00 . Witten argues that G exists because D has no zero modes and the 
operator ^ is formally (anti) self-adjoint (though see [40]). Any zero mode e of D must 
satisfy (B.2) with eoo = 0. Thus 


|Vef = 0, 


(B. 6 ) 


which implies e = 0 . 

Let us now consider the situation in the presence of the wall. As noted in section 2.2, 
the energy E then receives a contribution from the Brown-York stress tensor Tab [3] at the 
wall in addition to the familiar ADM-like boundary term E^ at infinity. We should thus 
impose a boundary condition at the wall such that this term appears in the analogue of 
(B.2). As noted in [41], this boundary condition should restrict half of the components of 
e. We therefore impose the condition 


P-e\ 


9So 


= 0 , 


(B.7) 
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where <9 Xq is the Dirichlet wall and 


p - = \{ 1 -• ( B - 8 ) 
This choice was also used in [41], It leads to 

(eoo4j£oo + [ T 0 AH A e > f |Ve| 2 > 0, (B.9) 

Jd So Jt, 

where the two terms on the left sum to the total energy (too Coc)E as desired. 

One would then like to show that there exists a Green’s function G = D —1 , where D is 
now the differential operator ^ restricted to spinors that vanish as r — > oo and satisfy (B.7). 
But the dT ,o boundary term in (B.9) is an obstacle, which in fact causes two problems. First, 
note that it arises from integrations by parts. This indicates that generally fails to be 
formally (anti-) self-adjoint under the boundary condition (B.7). This was not a problem 
in [41] where the analogue of dT,o was a surface on which Tab = 0. The dT ,o boundary 
term also obstructs the proof that V 7 has no zero modes, as instead of (B.6) one finds only 

/ T A h A £> [ |Ve| 2 >0. (B.10) 

Jd S 0 Jt. 

It thus appears that this method yields no positivity result for spacetimes with Dirichlet 
walls. Our original purpose in investigating this issue was to show that there is no tension 
between our negative kinetic energies and the E > Q argument of [29] which holds in 
spacetimes with no such walls. But the obstructions above also raise the question of whether 
the wall may in fact make the gravitational Hamiltonian unbounded below. 


C Convergence of the derivative of ip 


In the main text we showed that the sequence {Aip k } is uniformly convergent to zero as 
k -A oo. We now show that the same is true of derivatives d^Aip k and we place a bound 
on the rate of convergence, i.e. we derive (2.13). 

The proof of this convergence relies on the fact that is also a solution to the 

Laplace equation (since d commutes with the Laplacian) and the mean value property of 
harmonic functions. Using these two facts, and letting Atpk = ip — ipk, it follows that 


diAip k (x ) 


4nr 3 


' B(x,r) 


diAip k = 


Anr 3 


'dB{x,r) 


nAipk, 


(C.l) 


where B(x,r ) is a ball centered at x with radius r and r is unit normal to dB. Taking the 
absolute value of both sides, and using |fj| < 1 we obtain 


\diAip k (x)\ < -max(|AV>fc|) < a A (1 <5)(fe+1) , 

r dB \ r / 


(C.2) 


for k > 1, where the last inequality follows from (2.12). Iterating this procedure gives 


\B^Aip k (x)\ < 


max (|9 V1) A Mx)\) < a 


0"^- 5 )( fc +i), 


(C.3) 
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Now let y be an arbitrary point that is at least a distance h/2 away from the boundary 
dT,. Then, letting r = h/2 we obtain 


\d^AMy)\ 


(C.4) 


We will now use the analyticity of harmonic functions and (C.4) to obtain bound on 
\d^Ap k \ over all of E. Since d^ Aip k is a solution to the Laplace equation it is ana¬ 
lytic in the domain of convergence of it’s Taylor series, given schematically by 


9 (n) A Mv + e)=J2 


(e ■ d) m (d^Ap k (x)j 


m— 0 


ml 


(C.5) 


Using (C.4) we obtain 


\d (n) Ap k (y + e)\ < 


aA (i-*)(*+i) 
(h/6) n 


E 

_m =0 


1 [6\e 


ml V h 


(C.6) 


To cover all of E we only need to consider |e | < h/2 which gives 


0 n) A Mv + ?)\< 


aX (l-S)(k+l) 

(h/6) n 


E 

_m =0 


ml 


ae 3\( l S)(k+ 1 ) 
(h/6) n ' 


(C.7) 


Letting /3 := aer we obtain (2.13). 


D Constructing ip and K 

In this appendix we justify the use of K\. defined in (2.45), as an approximation of K. 
Examining (2.44) we see that we are only interested in the quantities d ■ K and \d x K |. As 
noted in the main text we can explicitly verify that ip\ := —d- K\ satisfies the boundary 
conditions (2.22) to order 0{a). So, it only remains to check that \d x K\\ faithfully 
approximates \d x K |. 

Without loss of generality we may write 


K = dU + d x C, 


(D.l) 


where d 2 U = — ip. Thus d ■ K = —ip, so all that remains is to construct a C such that 
d 2 I< = 0. 

A simple calculation reveals that 


$ x K = -d 2 {Cp) + 80- CP) = DP 
d 2 K = dp’ — d x (Dp), 


(D.2) 


where D is defined by the relation Dp = d 2 (Cp). So, we must now find a D such that 
d 2 I\ vanishes. We will show that such a D is given by 


B 2 (Dp) = 0, r-(—dp + dxDp) 




= 0, 


(D.3) 


26 - 











where r is the unit normal to 9X. To see that this is sufficient let L := d 2 K. We then have 

d ■ L = = 0, d x L = —^(D $) + d(d ■ D $) = 0, (D.4) 

where we have used d ■ K = — if} and d^D = 0. Therefore, L is the gradient of a solution to 
the Laplace equation, i.e. 

L = dV, &V = 0. (D.5) 


Furthermore, a complete set of boundary data for V is given by 


r ■ dV 




f • d 2 K 


8T, 

f ■ (—d'tjj + d x D B) 



(D.6) 


This boundary data is consistent with the solution V = 0 and so by the usual arguments, 
this must be the unique solution. Therefore 


iPK = 0 . 


(D.7) 


It is now straightforward to compute 

D 1 := 0 x K\ ) ■ B (D.8) 

and we find that D\ satisfies the boundary conditions associated with D up to terms of 
order 0(a). Therefore we obtain the result 

d-(K~ K\) ~ a, \d x (K - K\)\ ~ a. (D.9) 

It is straightforward to extend this construction to an arbitrary for which the errors 
described above would be 0(a k ), however higher order terms do not contribute to our final 
result. 
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